function res = FD_cos(omega,x,h)
res = (cos(omega*(x+h))-cos(omega*x))/h;